Origin of maximal symmetry breaking in even PT-symmetric lattices 
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By investigating a parity and time-reversal (VT) symmetric, A'^-site lattice with impurities ±27 
and hopping amplitudes to (it) for regions outside (between) the impurity locations, we probe the 
origin of maximal "PT-symmetry breaking that occurs when the impurities are nearest neighbors. 
Through a simple and exact derivation, we prove that the critical impurity strength is equal to the 
hopping amplitude between the impurities, 7c = tt, and the simultaneous emergence of A'^ complex 
eigenvalues is a robust feature of any "PT-symmetric hopping profile. Our results show that the 
threshold strength 7c can be widely tuned by a small change in the global profile of the lattice, and 
thus have experimental implications. 



Introduction: The discovery of "complex extension of 
quantum mechanics" by Bender and coworkers ^ [5] 
set in motion extensive mathematical i3j-t5| and theo- 
retical investigations [B] of non-Hermitian Hamiltonians 
H-pj- = K + V that are symmetric with respect to com- 
bined parity {V) and time-reversal (T) operations. Such 
continuum or lattice Hamiltonians [TlllOj usually con- 
sist of a Hermitian kinetic energy part, K — K\ and 
a non-Hermitian, T^T-symmetric potential part, V — 
VTVVT ^ Although it is not Hermitian H-pr 

has purely real eigenvalues E = E* over a range of pa- 
rameters, and its eigenfunctions are simultaneous eigen- 
functions of the combined T'T-operation; this range is 
defined as the T'T-symmetric region. The breaking of 
■pT-symmetry, along with the attendant non-reciprocal 
behavior, was recently observed in two coupled optical 
waveguides |111 [T^j and has ignited further interest in 
T'T-symmetric lattice models. These evanescently cou- 
pled waveguides provide an excellent realization [13j of an 
ideal, one-dimensional lattice with tunable hopping |14j . 
disorder |15| . and non-Hermitian, on-site, impurity po- 
tentials [THlfTT]. 

Recently nonuniform lattices with site-dependent hop- 
ping ta{k) — to [k{N — k)]"^^ and a pair of imaginary 
impurities ±«7 at positions (m, m) have been extensively 
explored [T7H^ . where fh ^ N + I ~ m and TV » 1 is 
the number of lattice sites. The 'PT-symmetric phase in 
such a lattice is robust when a > 0, the loss and gain 
impurities ±17 are closest to each other, and 7 < 7c 
where the critical impurity strength is proportional to 
the bandwidth of the clean lattice, 7c c>c Ato{N/2)" . For a 
generic impurity position m, when the impurity strength 
7 > 7c (ti) increases the number of complex eigenvalues 
increases sequentially from four to TV — 1 when N is odd 
and to N when it is even. In an exceptional contrast, 
when m = N/2 - nearest neighbor impurities on an even 
lattice - all eigenvalues simultaneously become complex 
at the onset of PT-symmetry breaking. This maximal 
symmetry breaking is accompanied by unique signatures 
in the time-evolution of a wavepacket |20j . 



lek@iupui.edu| 



These results raise the following questions: Is this ex- 
ceptional behavior limited to lattices with a-dependent 
hopping or is it generic? Which factors truly determine 
the critical impurity strength 7c(iV/2) in the exceptional 
case? How does the critical impurity strength 7c ('tt-) de- 
pend upon lattice parameters and impurity positions? 

In this Brief Report, we investigate an A^-site lattice 
with impurities ±27 at positions (m, m) and a constant 
hopping amplitude io('^b) for sites outside (between) the 
parity-symmetric impurity locations. Our two salient 
results are as follows: i) When m = N/2, we analyt- 
ically prove that all eigenvalues simultaneously become 
complex when 7 > 7c(A^/2) = tb- This robust result 
is true for any symmetric distribution of real hopping 
amplitudes, ii) When <{, 3> to, the critical impurity 
strength 7c(?7i) —J- ifc irrespective of the impurity posi- 
tion TO. When th < to, the critical impurity strength 
7c (m) ~ t^ where the exponent r]{d) ^ d increases mono- 
tonically with the distance c? = -|- 1 — 2to between the 
impurities. Thus, the P7~-symmetry breaking threshold 
can be substantially tuned without significant changes in 
the global hopping-amplitude profile of the lattice, and 
the exceptional nature of the m = A/2 case is due to the 
ability to partition the system into two, and exactly two, 
pieces. 

Tight-binding Model: We start with the Hamiltonian for 
a one-dimensional, tight-binding, non-uniform lattice 

AT-l 

(1) 

where ajj(a„) is the creation (annihilation) operator for a 
state localized at site n, and the hopping function is given 
by t{i) = t(, > for TO < i < TO — 1, and t{i) = to > Q 
otherwise. This Hamiltonian continuously extrapolates 
from that for a lattice of length d = A^ + 1 — 2to with 
impurities at its end when tf, ^ to, to that of a pair 
of disconnected lattices, one with the gain impurity and 
the other with the loss impurity, when t^ <^ to- Note 
that the critical impurity strengths in these two limits are 
known [171 EI] • Due to the constant hopping amplitude 
outside or between the impurity locations, an arbitrary 
eigenfunction \il>) = 'Yli^=i'4'{'^)'Aifi) with energy E can 
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be expressed using the Bethe ansatz as 

{Asin{kn), 1 < n < m, 

Psm{k'n) + Qcos{k'n), m, < n < fh, (2) 
Bs'm{kn), fh < n < N. 

Here E{k,k') — — 2tocos(/c) = — 2tbCos(fc') defines 
the relation between the quasimomenta k,k'. In the 
T^T-symmetric phase, the energy spectrum of Eq.Q is 
particle-hole symmetric [52], and the eigenenergies sat- 
isfy \E\ < 2max(io,ib)- Note that the relative phases 
of ipin) are the same at different points within each of 
the three regions, although there may be a phase differ- 
ence between wavefunctions in different regions. There- 
fore, without loss of generality, we may choose '4^{n) to 
be real for 1 < n < m. By considering the eigenvalue 
equation H-p-j-\ip) = E\ip) at points m,m + 1 and their 
reflection counterparts, it follows that the quasimomenta 
{k, k') obey the equation [3T] 

M{k, k') = [s\T? [k{m + 1)] + r2 siT?{km)\ 

X sin [k'{N + 1 - 2m)] + sm^{km) 
X sin [k'{N - 1 - 2m)] - 2TbSin{km) 
X sin[fc(m + l)]sin[fc'(Ar-2m)] =0, (3) 

where T — j/to and Tb ~ tt/to denote the dimensionless 
impurity strength and hopping amplitude respectively. 
Note that when 2mm{to,ti,) < \E\ < 2 max{to , tb) , k is 
real and k' is purely imaginary (or vice versa), whereas 
for \E\ < 2mm{to, tb), both k,k' are real. Thus, Eq.([3| 
represents two distinct equations in these two cases. 




FIG. 1. (color online) a) Left-hand panel shows dimension- 
less critical impurity strength Tc{d) = 7c/to as a function of 
dimensionless hopping amplitude <Tb = tb/to < 1 for vari- 
ous inter-impurity-distances d in TV = 20, 21 lattices; note the 
logarithmic scale. It follows that Tc{d) vanishes with a power- 
law behavior as Tb 0, as expected on physical grounds, b) 
Right-hand panel shows the critical impurity strength Tdd) 
as a function of Tj, > 1 for various values of d. Although at 
Tb = 1, the critical strength Tdd) reduces with distance d be- 
tween the impurities, for Tb > 2 the critical impurity strength 
Fc — >■ Tb (7c — >■ tb) irrespective of d. 



The right-hand panel in Fig. [T] shows the dimensionless 
critical impurity strength Tc{d) = 7c("^)/^o as a function 



of Tb = tb/tQ > 1 for various inter-impurity-distances 
d = N + 1 — 2m in an = 20 even lattice; we obtain 
similar results for an odd lattice. We find that 7c — > tb 
quickly for tb/to > 1; when tb/to ^ 1, the lattice reduces 
to one with d + I sites, impurities at its end points, and 
the result -fc = tb is expected [21]. The left-hand panel 
shows Tc{d) vs. Tb on a logarithmic scale in A*" = 20 
and iV = 21 lattices for Tb < 1. As the distance d be- 
tween the impurities increases, corresponding critical im- 
purity strength decreases as a power-law, T(.{d) cc T^^'*-* 
where the exponent r]{d) ^ d. This behavior can be 
qualitatively understood as follows: the system is in the 
T^T-symmetric region if the frequency ^ 7/^0 &t which 
particles are created at the gain-impurity site m is lower 
than rate at which these excess particles can hop over 
to the loss-impurity site, where they are absorbed at fre- 
quency J /to- Since tb is the hopping amplitude at 
sites between the impurities, it follows that the effective 
frequency of hopping from the gain- to the loss-site de- 
creases with d as Tl^. Indeed, when tb/to ^ 1, the system 
is divided into two, non T'T-symmetric, uniform lattices, 
one with the loss impurity and the other with the gain. 
It follows, then, that 7c — > as ff,/io — > 0. 
Origin of Maximal Symmetry Breaking: Now let us con- 
sider the m = N/2 case, where Eq.(|3]) reduces to 
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It follows from Eq.Q that the T'T-symmetry breaks 
maximally when 7 > 7c(A^/2) = tb and is accompanied 
by the simultaneous emergence of N complex (not purely 
imaginary) quasimomenta and eigenenergies. Since the 
bandwidth of the clean lattice is determined by both hop- 
pings (to, tb), it follows that the critical impurity strength 
is independent of the lattice bandwidth. 

To generalize this result, we consider the system with 
an arbitrary, PT-symmetric, position-dependent hop- 
ping profile ifc — tfj-k and real energy eigenvalues. 
Since the hopping and eigenvalues are real, the eigen- 
value difference equations imply that for any eigenfunc- 
tion we can choose the coefficients V'(^) to be real for 
1 < fc < m. A real eigenvalue e and the (real) coefficients 
a = (j){N/2) and /3 = (j){N/2 — 1) of its corresponding 
eigenfunction = J2iLi N) satisfy 



det 



tN/2-if3 + (e - 

tN/2Ct 



il)a tfq/2a 

tN/2-il3 + (e + «7)" 



0, 
. 

where we have used the T'T-symmetric nature of eigen- 
functions to deduce that 0(iV/2+l) e'^a, (j){N/2+2) = 
e*-^/3. Thus, when 7 > 7c = tjv/2 = tb, the eigenvalue e 
must become complex. Since this result is true for all 
eigenf unctions, it follows that the 'PT-symmetry breaks 
maximally and the critical impurity strength is solely 
determined by the hopping amplitude between the two 
impurities. This robust result also explains the fragile 
nature of T'T-symmetric phase in lattices with hopping 
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function ta{k) for a < [50]: in this case, the lattice 
bandwidth Aq ^ TV"'"'/^ whereas the hopping amph- 
tude between the two nearest-neighbor impurities scales 
as tb ^ A^~l"l. Therefore the critical impurity strength 
^c/ ^ A^^l"!/^ — > as — > oo. A similar analysis for 
closest impurities in an odd-A^ lattice shows that, due to 
the presence of a lattice site between the two impurity 
positions m = (A*" — l)/2 and in = {N + 3)/2, the cor- 
responding critical impurity strength 7c depends on the 



details of the eigenfunction. 

Thus, the maximal symmetry breaking only occurs in 
an even, 'PT-symmetric lattice with nearest-neighbor im- 
purities, and its origin is the ability to naturally partition 
such a lattice into exactly two components. 
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